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NOTATION 


Polynomial  coefficients 


Breadth  of  ship  section 


Breadth  of  flat  bottom 


Diameter  or  thickness  of  parallel  middle  body 


Diameter  or  thickness  of  flat  face 


Fullness  factor.  Equation  (26) 


Polynomial  associated  with  k 


Polynomial  associated  with  k. 


Curvature 


Rate  of  change  of  curvature  with  arc  length 


k at  x = 0 


k at  x = 1 


Polynomial  for  restraining  conditions 


Arc  length 


Axial  coordinate 


Axial  length  of  forebody 


x Normalized  axial  coordinate 

Y Offset  from  centerline 

y Normalized  offset 

Z Vertical  distance 

Vertical  height  of  curved  bilge 
a.  Adjustable  conditions 

g.  Restraining  conditions 
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ABSTRACT 


Polynomial  families  that  contain  a cube-root  term 
are  developed  to  satisfy  requirements  of  infinite 
slope  and  zero  curvature.  Adjustable  parameters  are 
used  to  provide  a wide  range  of  curves  without  In- 
flection points.  Applications  are  to  flat-faced 
underwater  bodies  and  wall-sided  ship  sections. 


ADMINISTRATIVE  INFORMATION 

This  work  was  authorized  and  funded  by  the  Naval  Sea  Systems  Command 
(Code  03512)  Task  Area  SR  023  01003,  Element  61153N. 

INTRODUCTION 

The  inability  to  obtain  infinite  slope  with  ordinary  polynomial 

expressions  has  led  to  addition  of  a square-root  term  in  describing  round- 
1 2 

nosed  bodies.  ’ When  zero  curvature  is  also  required,  the  square-root 
term  is  inapplicable.  It  is  now  proposed  to  substitute  a cube-root  term 
to  provide  both  infinite  slope  and  zero  curvature.  This  is  very  useful 
since  polynomial  expressions,  even  with  the  addition  of  a cube-root  term, 
lend  themselves  readily  to  analytical  manipulation.  Polynomials  remain 
polynomials  under  either  differentiation  or  integration. 

The  requirements  of  infinite  slope  and  zero  curvature  occur  for  under- 
water bodies  with  flat-faced  noses  and  for  curves  joining  ship  bottoms  to 
wall-sided  ship  sections,  among  others. 

The  requirements  of  infinite  slope  and  zero  curvature  have  also  been  met 

3 

by  previously  proposed  "cubic"  polynomials,  when  the  dependent  variable  is 
cubed  and  the  independent  variable  is  a polynomial  expression.  The  "cube 
root"  polynomial  now  being  proposed  may  be  more  convenient  to  the  designer 


Granville,  P.S.,  "Geometrical  Characteristics  of  Streamlined  Shapes," 
Journal  of  Ship  Research,  Vol.  13,  No.  4,  pp.  299-313  (Dec  1969).  A 
complete  listing  of  references  is  given  on  page  14. 

2 

Granville,  P.S.,  "Geometrical  Characteristics  of  Noses  and  Tails  for 
Parallel  Middle  Bodies,"  International  Shipbuilding  Progress,  Vol.  21, 

No.  233,  pp.  3-19  (Jan  1974). 

3 

Granville,  P.S.,  "Geometrical  Characteristics  of  Flat-Faced  Bodies  of 
Revolution,"  Journal  of  Hydronautics,  Vol.  7,  No.  4,  pp.  166-169  (Oct  1973). 
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since  the  dependent  variable  is  linear  instead  of  cubed.  Usage  will  deter- 
mine the  relative  merits  of  these  two  systems  in  hydrodynamic  applications. 

1-3 

It  has  been  found  that  two  adjustable  parameters  provide  a wide 
range  of  geometric  families  and  are  still  not  analytically  unwieldly  as  in 
determining  geometric  limitations.  Polynomial  expressions  for  one  adjust- 
able parameter  or  for  none  are  also  readily  obtainable  from  a two-parameter 
system. 

Without  loss  of  generality,  the  curves  to  be  considered  start  with 
infinite  slope  and  zero  curvature  and  end  with  zero  slope  and  zero  curva- 
ture. Other  conditions  may  also  be  considered  without  difficulty.  The 
two  adjustable  parameters  to  be  used  are  the  rates  of  change  of  curvature 
at  each  end  of  a curve. 

Normalized  coordinates  are  employed  and  the  resulting  polynomial  ex- 
pressions are  examined  for  suitable  ranges  of  values  for  the  two  adjustable 
parameters.  A particularly  stringent  geometric  requirement  is  absence  of 
inflection  points.  The  allowable  range  of  adjustable  parameters  is  ob- 
tained by  an  envelope  analysis. 

The  requirement  for  zero  curvature  at  both  ends  of  the  curve  is  used 
to  eliminate  discontinuities  in  curvature  at  junctions  with  curves  of  zero 
curvature.  Discontinuities  in  curvature  are  hydrodynamically  undesirable 
since  they  lead  to  pressure  changes  which  may  result  in  flow  separation 
and/or  cavitation. 


CUBE-ROOT  POLYNOMIALS 
In  general,  a cube-root  polynomial  y[x]  is 


1 n = N 


3 , Y1  n 
= a,  x + / ax 


n = 0 


where  a,  and  a are  coefficients  with  values  to  satisfy  boundary  conditions. 
1 n 


2 


Derivatives  are 


dy  _ 1 
dx  3 


(2) 


ll 

3 


x 


, i\  n-2 
n (n  - 1)  x 


(3) 


etc. 


At  x = 0, 


dx 


-*»  or  the  slope  is  infinite  as  required. 


Curvature  k may  be  expressed  as 


or 


(4) 


(5) 


It  follows  by  substitution  of  Equations  (2)  and  (3)  into  Equation  (5) 
that  k=0  at  x=0  as  required. 


TWO-PARAMETER  CUBE-ROOT  POLYNOMIALS 


The  convenient  adjustable  parameters  are  the  rates  of  change  of 
curvature  at  the  ends  of  the  curves.  Normalized  coordinates  0<x<l,  CKy£l 
are  to  be  used.  The  case  considered  here  will  have  zero  slope  and  zero 
curvature  at  x=l. 

From  the  curvature  given  by  Equatioi.  (5),  the  rate  of  change  of  curva- 

die  ^ 

ture  with  arc  length  s,  — = k is 


k 


_ dk 
ds 


(6) 


3 


dy  is  used. 


where  ds  = 


At  x = 0,  — 
dy 


0. 


Then 


(7) 


For  curvature  written  in  the  form  of  Equation  (4),  there  results 


_ dk  _ / d\* 


ds 


V 


dx“ 


-3 


-3 


(8) 


at  x = 1. 


dx  ,2 
dx 


0 


Then 


(9) 


In  the  case  of  being  considered,  the  cube-root  polynomial  has  two 
adjustable  parameters  0L,  i=l,  2,  the  rates  of  change  of  curvature  at  x=0 
and  x-l,  and  four  boundary  conditions  8 ^ , j=l,  2,  3,  4,  the  two  end  coordi- 
nates and  zero  slope  and  zero  curvature  at  x=l.  This  requires  a cube-root 
polynomial  of  the  fourth  degree 


y = 


a^x 


+ a + 


a^x 


2 . 
a^x  + 


3 . 4 

a~x  + a.x 
3 4 


(x0) 


4 


To  evaluate  the  coefficients,  the  adjustable  parameters  and  boundary 
conditions  Bj  are  substituted  into  the  cube-root  polynomial  equation, 
Equation  (10) 


(ID 


r 


ki  = 


4) 


or 


dx 


10 

27 


al  + 


6a^  + 


24a,  = k. 


x=l 


(12) 


B, : x = 0,  y = 0 or  a = 0 
1 o 


(13) 


B„:  x = 1,  y = 1 or  a,  + a + a,  + a_  + a.  + a = 1 
2 1 o 1 2 3 4 

3 

B3=  x = 1,  = 0 or  ~ a1  + ax  + 2a2  + 3a3  + 4a^  = 0 

3 

,2 

B^:  x = 1,  — ^ = 0 or  - -g  a^  + 2a2  + 6a3  + 123^  = 0 
dx  -r 


Solving  as  simultaneous  equations  results  in  a^ 

I 


(14) 


(15) 


(16) 


(17) 


Substitution  of  these  values  of 
in  the  adjustable  parameters  as 
the  form 


coefficients  and  gathering  of  terms  results 
coefficients  of  independent  polynomials  of 
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(18) 


'(I 


Kq[x]  + k1  Kl[x] 


Q [x] 


where 


Ko  = x3  + |-  (40x3  - 165x2  + 264x  - 220) 


(19) 


K,  = 


x , . , 3 

6 (x  “ X> 


(20) 


and 


Q = 1 - (x  - 1) 


(21) 


This  is  the  desired  two-parameter  cube-root  polynomial  satisfying  the 
specified  boundary  conditions. 


PERMISSIBLE  VALUES  OF  k AND  k, 

o 1 

Not  all  combined  values  of  kQ  and  k^  result  in  desired  shapes.  A 

stringent  limitation  is  to  have  no  inflection  points.  The  range  of  such 

values  is  given  by  a limiting  envelope  curve  which  is  derived  as  follows. 

,2  ,3 


Two  conditions  are  to  be  satisfied, 
two  simultaneous  equations 


_ 


= 0 and 


Ll.  = 


= 0,  which  result  in 


dx 


dx“ 


(r)  - 


| X 3 + |y  (80x2  - 165x  + 88) 


+ k^  (2x  - 3x  + 1) 


- 12  (x  - 1)  = 0 


(22) 


1 
3 r 


10 


|y  (x  + 32x  - 33 


+ k (4x  - 3)  -24 (x  - 1)  = 0 


(23) 


6 


(22x  - 24x  + 5)  + 3(30x  - 64x  + 33) 


where  x is  now  an  implicit  parameter.  For  the  range  of  values  0<  x _<  1, 
the  result  is  the  envelope  curve  of  Figure  1.  Values  of  the  adjustable 


6 \ 3 ^ 

parameters  [ — } and  k^,  which  give  noninflected  curves,  are  bounded  by 


the  axes  and  the  envelope  curve. 


RELATIVE  FULLNESS  OF  SHAPE 

A simple  measure  of  fullness  is  a prismatic  coefficient  F 


which  for  the  cube-root  polynomial,  Equation  (18),  becomes 


1 


o 


(26) 


(27) 


For  constant  values  of  F,  straight  lines  are  drawn  in  Figure  1.  Fuller 
shapes  are  in  the  lower  right-hand  part  of  the  figure. 


ONE-PARAMETER  CUBE-ROOT  POLYNOMIALS 


A one-parameter  family  of  cube-root  polynomials  may  also  be  developed 
1 

/ & \3 

if  ( — 1 , for  example,  is  designated  as  the  only  adjustable  parameter.  The 

\k  / 


result  is  a cube-root  polynomial,  one  degree  less  or  a cubic  cube-root 
polynomial.  In  this  case  is  zero,  and  Equation  (17)  gives 


o 


Figure  1 shows  the  equation.  It  is  to  be  noted  that  the  permissible  values 
of  and  for  the  two-parameter  system  are  much  greater  than  those  for 
the  one-parameter  system.  Substitution  of  Equation  (28)  into  the  two- 
parameter  polynomial,  Equation  (18) , produces  the  one  parameter  system 


y = \ z-J  ko[x]+q[x) 

^o 


where 


K = x3  - $(5x2  - 16x  + 20) 
o 9 


Q - 1 + (x  - 1)' 


d y 

A study  of  the  inflection  point  condition  — y = 0 indicates  a per- 

-v  dx 

missible  range  of  values  of  kQ  given  by 


21  </6_V  < 81 

16  “At  ) 

o 


for  no  inflection  points.  Figure  2 shows  the  curves  with  these  limiting 
values . 
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ZERO  SLOPE 


ZERO  CURVATURE 


0.2  0.4  0.6 


Figure  2 - One-  and  Zero-Parameter 
Cube-Root  Polynomials 

(Normalized  Coordinates) 


ZERO-PARAMETER  CUBE-ROOT  POLYNOMIAL 
A nonadjustable  cube-root  polynomial  may  also  be  designated  where  kQ 
and  are  not  adjustable  parameters.  In  this  case  a^  and  a^  are  zero, 

1 

/ 6 V 9 ~ 2 

and  from  Equation  (17)  it  follows  that  I — I = — and  k^  = — ; this  point  is 

^o 

plotted  in  Figure  1.  The  cube-root  polynomial  then  reduces  to 


9 3 . x 

y=-x  - x + — 


which  does  not  have  an  inflection  point  for  x < 1.  A plot  of  the  zero- 
parameter  cube-root  polynomial  is  shown  in  Figure  2. 


FLAT-FACED  NOSES 


> ! 


The  cube-root  polynomial  developed  here  may  be  used  to  bridge  a flat- 
faced nose  and  a parallel  middle  body  as  shown  in  Figure  3.  Both  two- 
dimensional  and  axisymmetric  bodies  may  be  accommodated.  If  is  the 


width  or  diameter  of  the  flat  face,  and  D is  the  width  or  diameter  of  the 
parallel  middle  body,  the  normalized  coordinates  are 


x = 


(34) 


and 


2Y  - D, 


y = 


D-D, 


(35) 


where  X is  longitudinal  or  axial  distance  from  the  flat  face 


is  length  of  the  nose 


Y is  normal  or  radial  distance  from  the  centerline 


The  origin  of  the  bridging  curve  is  at  the  outer  edge  of  the  flat  face, 
and  the  terminal  is  the  junction  with  the  parallel  middle  body.  At  both 
junctions,  zero  curvature  is  satisfied. 


WALL-SIDED  SHIP  SECTIONS 


I 


The  cube-root  polynomials  developed  here  may  also  be  used  to  bridge 
wall  sides  and  flat  bottoms  of  ship  sections  as  shown  in  Figure  4.  Here 
the  normalized  coordinates  are  given  by 


and 


B - 2Y 


(36) 


(37) 


where  Y is  the  offset  from  the  center  plane 
B is  the  breadth  of  the  ship  section 
B^  is  the  breadth  of  the  flat  bottom 
Z is  the  vertical  distance 
Z^  is  the  vertical  height  of  the  bridging  curve 

The  origin  is  the  junction  with  the  wall  side,  and  the  terminal  is  the 
junction  with  the  bottom.  At  both  junctions,  zero  curvature  is  satisfied. 


I 


Previous  efforts  to  use  ordinary  polynomials  to  fit  wall  sides 
resulted  in  approximate  expressions  with  terms  to  the  200th  power.  The 
cube-root  polynomial  serves  the  same  purpose  quite  readily  in  an  exact 
fashion. 


^Kerwin,  J.E.,  "Polynomial  Surface  Representation  of  Arbitrary  Ship 
Forms,"  Journal  of  Ship  Research,  Vol.  A,  No.  1,  pp.  12-21  (Jun  1960). 
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